Superconductivity and Field-Induced Magnetism in Pr^^Ce^CuC^ Single Crystals 
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We report muon-spin rotation/relaxation (/iSR) measurements on single crystals of the electron- 
doped high-T c superconductor P^-zCe^CuCU. In zero external magnetic field, superconductivity is 
found to coexist with Cu spins that are static on the ^tSR time scale. In an applied field, we observe 
a /i + -Knight shift that is primarily due to the magnetic moment induced on the Pr ions. Below the 
superconducting transition temperature T c , an additional source of static magnetic order appears 
throughout the sample. This finding is consistent with antiferromagnetic (AFM) ordering of the 
Cu spins in the presence of vortices. We also find that the temperature dependence of the in-plane 
magnetic penetration depth X a t in the vortex state resembles that of the hole-doped cuprates at 
temperatures above ~ 0.2 T c . 

PACS numbers: 74.25.Ha, 74.72.-h, 76.75.-M 



While there now exists a large body of convincing ex- 
perimental work on high-T c superconductors with hole- 
type carriers, the intrinsic properties of the electron- 
doped cuprates -rt^-zCe^CuC^ (R = La, Pr, Nd, Sm 
or Eu) have remained very elusive. Part of the problem 
stems from the difficulty of preparing single phase super- 
conducting samples. An additional complication is the 
presence of both Cu and rare-earth (i?) moments. The in- 
terplay between these two magnetic sublattices has been 
studied extensively in the undoped parent compounds 
.R2CUO4 EI- In Ce-doped samples, there is the possibil- 
ity that the magnetic exchange interactions involving the 
Cu02 planes play an important role in superconductiv- 
ity. There is also some uncertainty about the symmetry 
of the superconducting order parameter. Although, re- 
cent phase sensitive Q and angle-resolved photoemission 
experiments are consistent with d-wave pairing sym- 
metry, measurements of thermodynamic quantities such 
as the in-plane magnetic penetration depth A a b, generally 
do not show the expected linear temperature dependence 
at low T that is characteristic of a "clean" superconduc- 
tor. Many of the early results for A a &(T) in electron- 
doped cuprates were obtained in the Meissner phase of 
Nd2- a; Ce 2 ;Cu04 (NCCO) thin films and single crystals. 
While these results favored s-wave symmetry 0, 0, @j 
it has since been recognized that the paramagnetism of 
the Nd ions strongly affects the low-T behavior of \ a b(T) 
000 El ■ Consequently, attention has shifted primar- 
ily to the Pr2-ajCe a; Cu04 (PCCO) system, where crystal 
electric field (CEF) splitting of the Pr 3+ J = 4 manifold 
results in a nonmagnetic singlet ground state. Unfortu- 
nately, more recent measurements of \ a b(T) in PCCO 
thin films and crystals using a variety of techniques have 
not yielded consistent results @, H EE IH 13 111 El ■ 



To avoid extrinsic effects related to the sample surface, 
we have utilized the /iSR technique. /iSR is a local probe 
of the bulk and is one of the few techniques that can 
measure X a b(T) in the vortex state — sensitive to the su- 
percurrents circulating around the vortex cores, rather 
than the screening currents which flow near the sam- 
ple surface. Attempts to accurately measure X a b(T) in 
NCCO by /iSR have been precluded by the considerable 
Nd-moment contribution to the muon-spin depolariza- 
tion rate EH- In this Letter, we report on the first /LiSR 
study of superconducting PCCO single crystals. 

We studied three single crystals of PCCO grown by 
a directional solidification technique in AI2O3 crucibles 
using a CuO-based flux [lq. Reduction of oxygen to 
induce superconductivity was achieved by encapsulating 
each single crystal in polycrystallinc PCCO in the pres- 
ence of flowing Ar at 900-1000 °C (Ref. 17]). Here we 
report on representative data from one of the single crys- 
tals, noting that qualitatively similar results were ob- 
tained in a partial /iSR study of the other two. The 
crystal was 0.07 mm thick, having an a-b plane surface 
area of ~ 6.5 mm 2 and a mass of 3.74 mg. Although 
resistivity measurements indicate a T c value of 25(1) K, 
bulk susceptibility x" measurements show that the dia- 
magnetic signal assumes a constant value below ~ 16 K 
(see Fig. 1 inset). 

The measurements were carried out on the M15 and 
M20B positive muon beam lines at TRIUMF, 

Canada. In a /iSR experiment the spin of the implanted 
fi + precesses at a frequency uj^ = 7 m S m , where is the 
local field at the /i + site and 7^ = 13.5534 kHz G _1 is 
the muon gyromagnetic ratio. The PCCO single crystals 
are the smallest samples studied to date by conventional 
fj.SH. Consequently, a special arrangement of scintilla- 
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FIG. 1: Time evolution of the asymmetry (i.e. A(t) = AP(t), 
where A is the signal amplitude) in zero external field at T= 
80 K, 40 K and 2.3 K. The initial muon-spin polarization was 
perpendicular to the c-axis. Inset: Temperature dependence 
of the resistivity (solid diamonds) , and the bulk susceptibility 
at H II c = 1 Oe under FC and ZFC conditions. 
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FIG. 2: FFTs of the muon-spin precession signal in PCCO. 
Below T c , the measurements were taken in an arbitrary se- 
quence as a function of T under FC conditions. The dashed 
vertical line indicates the value of the external field H = 91 Oe. 
Inset: Asymmetry spectra taken at T = 2.3 K under both FC 
and ZFC conditions. 



tor detectors inside a He gas-flow cryostat was needed to 
eliminate the contribution from muons that did not stop 
in the sample. Unfortunately, this setup is not possible in 
a dilution refrigerator, so 2.3 K is the lowest temperature 
that could be reached. Measurements were taken both 
in zero external field (ZF) and in a transverse field (TF) 
geometry with H directed perpendicular to the Cu02 
planes. 

In ZF, the muon-spin polarization Pit) is relaxed by 
magnetic moments of electronic and nuclear origin. As 
shown in Fig. 1, the relaxation rate of the ZF-/iSR sig- 
nal in PCCO increases with decreasing T. The absence 
of discrete frequency oscillations indicates that there is 
no onset of spontaneous magnetic order at temperatures 
above 2.3 K. Additional measurements carried out in a 
longitudinal-field geometry indicate that the muon spin 
is depolarized by randomly oriented internal fields which 
are static on the /j,SR time scale. Since the Pr ions have 
a nonmagnetic singlet ground state and the nuclear con- 
tribution to Pit) is small, we attribute the source of the 
magnetism to disordered Cu spins. 

TF-/iSR and magnetization measurements were car- 
ried out under both field-cooled (FC) and zero-field 
cooled (ZFC) conditions. Typical asymmetry spectra are 
shown in the inset of Fig. 2. In the ZFC case, pinning at 
the sample edges prevents flux from entering the bulk at 
T = 2.3 K and H < 300 Oe. Thus, the ZFC time spectrum 
in Fig. 2 resembles that observed in ZF. We note that 
this is also the case in high-quality YB^CuaOy-^ single 
crystals. Above 300 Oe, flux fully penetrates the sample 
at T = 2.3 K, but pinning centers prevent the formation 
of an equilibrium vortex lattice configuration. On the 
other hand, a well-ordered vortex lattice is achieved in 
the FC procedure. Fast Fourier transforms (FFT) of the 
corresponding muon-spin precession signals at H = 91 Oe 
are shown in Fig. 2. The FFT provides an approximate 



pict ure of the internal magnetic field distribution n(B) 
|18|. Above T c , the linewidth is due to the electronic and 
nuclear magnetic moments. Below T c , the lineshape is 
further broadened and becomes asymmetric as a result 
of the inhomogeneous field distribution created by a lat- 
tice of vortices. Because the contribution of the magnetic 
moments to the width of n(B) increases with increasing 
if, we have restricted the current study to low magnetic 
fields. 

The stopping site of the [i + has previously been identi- 
fied in the related undoped parent compound Nd2Cu04 
as near an O atom midway between adjacent Cu02 layers 
[T^| . Consistent with this site assignment, the TF-//SR 
signal in PCCO shows a single well resolved signal with 
an average frequency shifted relative to the Larmor pre- 
cession frequency of fi + in vacuum. The corresponding 
/i + -Knight shift is 

where Bq is the average internal field sensed by the 
muons, the second term is a correction for Lorentz and 
bulk demagnetization fields, p mo \ = 0.01744 mol/cm 3 is 
the molar density of PCCO, and N ks 1 is the demag- 
netization factor for a thin platc-like crystal. Figure 3 
shows that scales linearly with x" above T c , indi- 
cating that the relative frequency shift is not induced 
by the /i + . In the normal state there are several con- 
tributions to ifji. Susceptibility measurements have es- 
tablished that the field-induced moment on the Pr ions 
dominates the magnetic response of PCCO, which de- 
pends little on Ce-doping [2(| . Furthermore, as shown in 
the inset of Fig. 3, x is highly anisotropic due to strong 
crystalline electric-field (CEF) effects. For H || c, there 
are smaller contributions to from the tilting of the 
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FIG. 3: Kf, vs. x" at H = 3 kOe. The solid line is a fit 
to Eq. (|5J. Left inset: temperature dependence of the bulk 
susceptibility at H — 10 kOe applied parallel (x") and per- 
pendicular (x ± ) to the c-axis. Right inset: magnetic field 
dependence of Bo — fioH at T = 5 K and 25 K. Each data 
point at T=5 K represents a single FC measurement. 



Cu moments out of the a-b plane, Cu 2+ -0 2 ~-Pr 3+ su- 
pcrexchangc interactions and from the polarization of the 
conduction electrons |2^. The Pauli paramagnetic sus- 
ceptibility xo of the conduction electrons gives rise to 
a T-independent Knight shift, Kq = Aqxo, where Aq is 
the contact hyperfine coupling constant. The field in- 
duced moment on the Pr ions results in a T-dependent 
Knight shift Kj consisting of two contributions: (i) the 
dipole-dipole interaction between the localized Pr 4/- 
moments and the fj, + , and (ii) an indirect Rudermann- 
Kittel-Kasuya-Yosida (RKKY) interaction, producing a 
spin polarization of the conduction electrons at the /x + 
site. Assuming the Cu-moment contributions to be neg- 
ligible, ifj] at the axial symmetric /i + -site is given by 



o + (A c + A|f p ) x } ; 



(2) 



where A r and A 



dip are the contact hyperfine and dipolar 



coupling constants, respectively, and x^ ~ X ~ Xo- At 
the /i + site we calculate A%? =— 663.6 Oe/^B- The solid 
line in Fig. 3 is a fit to Eq. ©, yielding A c = 3231 Oe//z B 
and K > -4480 ppm. 

Below T c , there is a substantial increase of Bo that 
nearly coincides with the increase of the diamagnetic sig- 
nal under FC conditions (see Fig. 4). This cannot be ex- 
plained by the reduction of xo that arises from the forma- 
tion of Cooper pairs. We first considered the possibility 
that the increased field is a manifestation of the so-called 
paramagnetic Meissner effect |22^ . but the diamagnetic 
response observed in both the FC and ZFC magnetiza- 
tion does not support this interpretation. A second pos- 
sible origin, thought to be responsible for magnetorestric- 
tion enhancements [2^ and large diamagnetic shifts of B 
[24| in the superconducting state of the iJI^CuaCv-^ 
system, is the i?-moments induce screening currents in 
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FIG. 4: B (T)-B$ vs. B$ at T = 5 K (solids circles) and 
10 K (open circles). The solid curves are described in the text. 
Inset: Temperature dependence of Bo and x" taken under FC 
conditions. 



the Cu02 layers. However, here the Pr-moment is in- 
duced by the applied field, and thus is reduced below 
the normal-state value in regions well beyond the vor- 
tex cores. This is inconsistent with the FFTs in Fig. 2, 
which imply that all muons stopping in the sample see 
an increase in the local field. A third possible origin is 
that the vortices nucleate static Cu-spin order, as has 
been observed in the superconducting state of under- 
doped La2~ 2: Sr :E Cu04 [23. We note that while the dif- 
ference between Bq in the superconducting and normal 
states is reduced with increasing H, this does not im- 
ply that the additional source of field below T c weakens. 
For example, the solid curve in Fig. 4 is a fit assuming 
B (T) = [{B^) 2 + {B 1 :) 2 ] 1 I 2 , where B^ is the average in- 
ternal field in the normal state and B± «43 G and 28 G, 
at 5 K and 10 K, respectively. This indicates that the 
additional field B± that appears at the ^j + -site is pri- 
marily directed in the a-b plane. Calculations show that 
this is consistent with the La2Ni04 spin-structure that 
has been observed with /^SR in the parent compound 
Nd2Cu04 at low T 19]. Thus, long-range order of the 
Cu spins appears to be stabilized in PCCO by a low den- 
sity of vortices, perhaps due to the close proximity of the 
superconducting and AFM phases. 

The in-plane magnetic penetration depth A a {, can be 
determined from the TF-/1SR time spectra as described 
in Ref. 01- I n the normal state, the muon-spin pre- 
cession signal is well described by the polarization func- 
tion P(t) = e~( CTt ) cos(a; M £ + 0), where (3 w 1.6 and a 
is the temperature-dependent depolarization rate. The 
TF-/J.SR time spectra below T c were well fit to 

P(t) = e-^t) I nv (B) cos( 7MJ Bi + cj))dB , (3) 



where n v (B) is the field distribution of an hexagonal vor- 
tex lattice derived from Ginzburg-Landau theory, and 
er v is the depolarization rate due primarily to the effects 
of the magnetic moments — which occur on a length 
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FIG. 5: Normalized magnetic penetration depth A„ b (0)A^ b (T) 
as a function of T/T c for single crystals of PCCO at H = 90 Oe 
(sol id circles), I^-iSr^CuCU at H = 2 kOe (open squares) 
filH and YBa 2 Cu 3 06.95 at H=5 kOe (open circles) |2g. The 
solid curve is a guide for the eye. Inset: A^ b 2 vs. T in PCCO. 
The solid curve is a fit described in the text. 



scale that is 10 3 times smaller than A a &. Fits at low T 
which explicitly included a variable transverse field com- 
ponent B±, gave fi ~ 1.2 and a diamagnetic shift for 
the c-axis component of field. However, \ a b(T) could 
not be reliably obtained in this manner due to the de- 
creased sensitivity to B± with increasing T and the com- 
petition with other fitting parameters. Our analysis us- 
ing Eq. © shows that a v scales as 1/(T — 0) , where 
9 = — 18(5) K. At low H, only a small fraction of the im- 
planted muons stop near the vortex cores. As a result, 
the analysis is not very sensitive to the value of the coher- 



ence length £ a ;,. For example, increasing £ a fc from 30 A to 
60 A changes A a b(0) by less than 2 %. Figure 5 shows 
the temperature dependence of A~ b 2 . The absence of data 
below T — 2.3 K forbids an accurate determination of 
the limiting temperature dependence of A~ b 2 (T). How- 
ever, above this A~ b 2 (0)/A~ b 2 (T) shows reasonable agree- 
ment with that determined by /iSR in the hole-doped 
systems I^-zSr^CuC^ and YBa2Cu30y_5. The solid 
curve in the inset of Fig. 5 is a fit of the PCCO data to 
A; 6 2 (T) = A- 6 2 (0)[1-(T/T C ) 2 ], yielding T c = 15.9(2) K and 
A ab (0) = 3369(73) A. The values of T c and A ab (0) indicate 
that the bulk of our sample is primarily underdoped. For 
example, Meissner state measurements in Ref. | 13| on an 
underdoped PCCO film reported A a6 = 3100 A. This same 
study showed evidence for a crossover from s-wave to d- 
wave behavior in the underdoped regime. 

In conclusion, we have observed the onset of static 
magnetic order in the superconducting state of under- 
doped PCCO single crystals with a field applied perpen- 
dicular to the Cu02 plane. This result is consistent with 
vortex-induced AFM ordering of the Cu spins through- 
out the sample. In addition, our measurements of X~ b (T) 
in the vortex state appear consistent with hole-doped 
cuprates down to ~ 0.2 T c . However, a unique determi- 
nation of the low-temperature behavior by /iSR awaits 
the availability of larger single crystals. 
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